Thermal entanglement in a molecular chain 
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We investigate entanglement in a linear chain of A'' polar molecules coupled by dipole interaction. 
In our model, nearest neighbour interaction predominate, and we compute entanglement with the 
help of a two-party correlation entanglement measure. We find that, in this system, only excited 
states are entangled. Moreover, when an electrostatic field is applied, energy levels crossings occur, 
changing significantly the system's entanglement properties. We make a systematic study of the 
entanglement dependency on the inter molecular distance separating pairs of molecules, different 
partitions of the chain and physical parameters as the temperature and the electrostatic field's 
intensity, showing that it persists for relatively high temperatures and changes its nature with 
varying field. 
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The physics of quantum many-body interacting sys- 
tems is a passionating research field in the frontier of 
quantum mechanics, statistical physics and quantum 
field theory. About 30 years ago, quantum many-body 
systems have also been implicated in the domain of quan- 
tum information theory. At its origin itself quantum 
computers were conceived as physical systems allowing 
for the simulation of complex quantum many body sys- 
tems in a controllable context, working as analogue com- 
puters. Since then, the interlace between quantum many- 
body systems and quantum information, has increased 
and showed itself as a new manner of looking at quantum 
phase transitions (QPTs)[2|, as well as a major applica- 
tion of controllable physical systems, as trapped ultra- 
cold atoms and molecules |j] . Pioneer works studied 
the behaviour of a particular entanglement measure, the 
concurrence C Q, close to a QPT in a chain of spins 1/2 
interacting via the Ising Hamiltonian showing that 
singularities in the derivative of C also appear in the ther- 
modynamic limit. Long range spin models were equally 
studied from the entanglement point of view 0], and 
it was shown that discontinuities also occur for entan- 
glement measures close to QPT points. An application 
of such studies concerns collective physical properties of 
spin systems, as magnetic susceptibility, that can be di- 
rectly influenced by the presence of entanglement Q . 
In parallel to that, an important step has been taken in 
the domain of simulation of quantum many-body mod- 
els with the advent o f op tically trapped cold atom and 
molecular physics 0, llOl . HH . Theoretical models show 
that ultra cold atoms posed in optical lattices can simu- 
late Bose-Hubard spin models Experimentally, the 
observation of the Mott-insulator transition marked the 
beginnings of the domain of solid state physics with cold 
atoms [J • Technical developments turned possible 
the cooling and trapping of more complex structures, as 
diatomic molecules |13| . that present some potential ad- 
vantages with respect to atomic systems: polar molecules 
can have a strong dipole moment, controlled by manipu- 



lating the molecular internal states. This can be used as 
a toolbox to simulate sohd state spin models 

In the present paper, we develop a new way of looking 
at trapped cold polar molecules, different from the quan- 
tum simulation perspective. Since an ensemble of polar 
molecules is naturally a many-body interacting system, 
we describe its Hamiltonian's properties using techniques 
similar to those employed for treating entanglement in 
solid state spin systems. We deal here with a high spin- 
like system: the value of j, the single molecule angu- 
lar momentum, determines the dimension of each "spin" . 
Our system interacts at its own manner, via dipole-dipole 
coupling, giving rise to a qualitatively different behaviour 
than the one observed in condensed matter systems with 
respect to its singularities. In the present paper, we are 
mainly interested on the entanglement in different types 
of eigenstates of the molecular ensemble and how it can 
be modified by physical parameters as the intensity of an 
applied electrostatic field or the temperature T. We also 
analyse how the variance of the z projection of the total 
angular momentum J^, (AJ^)^, is affected by the same 
parameters. This quantity is analogous to the variance 
of the magnetisation in spin systems, and can be used to 
label different types of entangled states, as seen in the 
following. Our results have applications in the physics of 
cold molecules, helping the preparation and the charac- 
terisation of different collective states. From the funda- 
mental point of view, it exhibits a number of interesting 
features, as the presence of thermal entanglement 14| . 
long range correlations and a discontinuous behaviour. 

We deal here with a one dimensional chain of diatomic 
polar molecules coupled by dipole interaction. A simpli- 
fying assumption made here is that molecules arc suffi- 
ciently far apart, so that the dipole interaction can be 
treated as a perturbation coupling nearest neighbouring 
molecules only. The dipole interaction scales as 
where r is the inter-molecular distance, considered as 
constant. This means that interaction strengths lower 
than 1 /8th of the one coupling nearest neighbours are dis- 
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regarded. As an example, we can take the KRb molecule, 
with n = 1.2D and B w lOGHz. When, in a chain, 
KRb molecules are spaced of r « 5nm, we have that 
Vd,p/2B « 5%. 

We now move to the Hamiltonian of the system, which 
determines all the properties we wish to investigate. Each 
molecule is treated as a rigid rotor and the i-th molecule 
rotational eigenstate is |ji,mi), where ji is the total 
angular momentum quantum number and nii its pro- 
jection onto the z axis, taken as the chain axis. The 
free molecule's energy is Ei = + 1)B, where B is 
the rotational constant. The total energy of a chain of 
N free molecules in the eigenstate |ji,TOi) (2) ^2,^2) (2) 
... (g) \jN-i,mN-i) <E) \jN,mN) is E ^ BY,f + 1). 
We have here a highly degenerated system: for each 
molecule, for a given ji, there are 2ji + 1 states with 
different m^. When considering a chain of N molecules, 
we add degeneracies coming from all the permutations of 
the different jiS with the same energy. 

We now describe the dipole interaction, that may 
lift part of the degeneracy. By supposing that 
all the molecules have the same dipole moment /2, 
the dipole-dipole interaction between two neighbouring 
molecules is Vdip = |/xp/(27reor^)(— 2 cos6'i cos6'i+i + 
sin^i sin^i4_i cos (0i — ^i+i)), where 9i is the angle the 
i-ih molecule makes with the z axis and ipi is the 
corresponding azimuthal angle. Vdip possess sym- 
metries and selection rules simplifying the treatment 
of the problem. For a given pair of molecules, 
the dipole interaction has a non null matrix element 
{f^,m[\{j[J^^,m[_^_^\Vd^p\i^,m,)\j^+l,m^+l) only if = 
ji ± 1 and rrii + rrii+i = m[ + ni[_^-^. The inter- 
est of Vdip is that it creates entanglement. To see 
how it happens, we consider the simple case of two 
molecules, that helps giving an insight of what happens 
for ^ 2 molecules. The free Hamiltonian ground 
state is trivial, non-degenerated, and non-interacting. 
We then discuss the first excited state, 6 times degen- 
erated, of energy E = 2B, and composed by states 
|1,0)|0,0), |1,1)|0,0) and |1,-1)|0,0) and aU the other 
three permutations. Vdip creates a new eigensystem, 
composed of states |*^) = l/^/2(|l, 0) |0, 0) T |0, 0) |1, 0)), 
of energy f j."' = 2B ± \^\^ / {Gneor'^) and the degener- 
ated states 1*^) = 1/V2(|1,1)|0,0) ± |0;,0)|1,1)) and 
|*±i) = 1/a/2(|1,-1)|0,0) ± |0,0)|1,-1)), of energy 
= 2B± |^|V(127reor3). The energy difference be- 
tween both types of states occurs because for l^fj.) states, 
the angular probability distribution is maximum on the 
z axis, corresponding to 9i = 0, minimising Vdip. In 
contrast, for states, the maximum lies in a plane 

orthogonal to the z axis. The first excited state of 
the dipole interacting system, |^'° ), is pure and entan- 
gled. To measure the amount of entanglement of this 
state, we use the logarithmic negativity [l5|, defined as 
£ = log2 {2Af + 1), where J\f is the negativity, the sum of 



the negative eigenvalues of the partial transposed density 
matrix of the two party system [I^. In such effective 4 
level system, £ 7^ is a sufficient condition for detecting 
entanglement. We find that C = 1 for We can 

also calculate (AJ^)^ for this state, finding (AJ^)^ = 0: 
the angular momentum is orthogonal to the z direction, 
as mentioned previously. The next lowest energy states 
arc doubly degenerated: 1^'^"'^), with energy To 
study their physical properties we assume that they form 
an equally weighted mixed state. The resulting state, 
even if mixed, is also entangled, with C = 0.77. We also 
have {AJz)^ = 1, essentially different from the previous 
case. We can thus distinguish the two lowest energy one- 
excitation states from the point of view of entanglement, 
(AJz)^ and purity, showing that two classes of states can 
be defined in the subspace with one rotational excitation. 

We now discuss qualitatively the case of a chain of N 
molecules before moving to our results. Initially, the sys- 
tem's degeneracy of the first excited level is 3A^. Among 
this states, we distinguish A^ with one excitation of the 
type 1 1,0) (called from now on the Ho subspace) and 
2A with one excitation of the type |1,±1) (called from 
now on the Hi subspace). When taking Vdip into ac- 
count, states in the Ho subspace become completely non- 
degenerated, forming A^ different entangled states. In 
contrast, since Vdip is independent on the sign of m^, each 
energy level in the Hi subspace is two times degenerated, 
so that entanglement in this subspace is essentially dif- 
ferent than in Hq. Also, (AJz)^ = 1 always in Hi and 
null in 7^0- 

From this basis, we study now how the spectrum of 
the system, as well as its properties, can be modified 
by the application of an electrostatic, linearly polarised, 
uniform field in the z direction, interacting independently 
with each molecule. The field's interaction Hamiltonian 
is Ve ~ — ^ii?z cos^i, where E^ is the field ampli- 
tude. We start by analysing its effects over each molecule 
independently, ignoring Vdip- Due to the field's selec- 
tion rules, single molecule states satisfying j,- ~ ji ± 1 
and m'i = mi only are coupled. However, instead 
of considering the exact field's eigenstates, we project 
them into the subspace that is coupled at first order by 
Vdip- In this subspace, the field's eigenstates are |+)j = 
cosi/3i|l,0)i-|-sin y.t|0,0), \-)i = cos y^lO, 0)t -sin ipt\l, 0) , 
with cosifii = y/{B + X)/{2\), and A = ^/W^+ElJi- 
The total energy of |-|-) and |— ) is E± = B ± X- States 
of the type |1,±1) remain unchanged (after projection) 
but their energy is shifted. Such approximation has been 
corroborated by including the ji = 2 states in numerical 
calculations. 

In a second step, we consider a molecular chain, and 
compute the effects of Vdip- To do so, the best strat- 
egy is to express Vdip in the basis of the field eigen- 
states. In this case, the system's Hamiltonian is tri- 
diagonal and can be analytically solved We find 

then that the ground state is still non entangled in first 
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order. The first excited level has TV — 1 molecules in 
state |— ) and one molecule in one of the three possible 
excited states, |+), |1,±1). Again, two subspaces can be 
defined, the one with one |+) excitation (7i+) and the 
one with one |1, ±1) excitation {Hi). The system's spec- 
trum for = 50 depends on the dimensionless electric 
field amplitude = E z 4:^/371 Cor^ /In] as shown in Fig. 
[TJ There are no level crossings inside a same subspace, 
since states in a same subspace are coupled to the field 
in the same way. On the other hand, the Tii subspace is 
differently coupled to the field than the H+ one, leading 
to level crossings between the two subspaces. Because 
each subspace has different values, experimental char- 
acterisation of each subspace can be done by polarisation 
dependent spectroscopy. 

We study now the observable and relevant conse- 
quences such level crossings have in the system's entan- 
glement. We use here two methods to calculate entan- 
glement. In the first one, the entanglement of a pair of 
molecules is evaluated by computing £ after tracing out 
the other N — 2 molecules of the chain. We can then 
obtain Ld, the total value of £, summed over all pairs 
of molecules at a same relative distance d. This is use- 
ful for giving information on the range of entanglement 
and on how it is distributed over the chain. A second 
method to compute entanglement consists of calculating 
L'p, the entanglement of one molecule to the rest of the 
chain. This quantity, in this finite model, depends on 
the chosen molecule's position p: it is maximised for the 
central molecule and minimised for the extrema ones. We 
observe that states in the same subspace share the same 
value of entanglement using such definitions. In Fig. [2] 
we show Ld as a function of Cz for the lowest energy ex- 
cited level in a chain of A' = 50 molecules. We see that 
at the field amplitude corresponding to the energy cross- 
ing point, entanglement suffers a discontinuous change. 
Also, for each class of state, entanglement does not sig- 
nificantly depends on e^. We observe that entanglement 
decreases as the distance between pairs of molecules d 
increases, becoming close to zero at the extrema of the 
chain. Discontinuity in is observed for all values of 
d, which is a particularity of this system when compared 
to other nearest neighbours models Q. In the same fig- 
ure, Ljg and L[ are plotted for a chain of 50 molecules, 
and both quantities suffer a discontinuous change at the 
energy levels crossing point. 

In the present system, entanglement is different from 
zero only for excited states. Consequently, it displays 
thermal entanglement: starting from the ground state, 
entanglement is created by raising the temperature T 



14j . A natural question is how entanglement in thermal 



equilibrium depends on and on a re-scaled tempera- 
ture T = UbT / B , where ks is the Boltzmann's constant. 
For low enough values of T, it is reasonable to suppose 
that only the subspace with one rotational excitation is 
populated. Fig. [3] shows as a function of T and 




FIG. 1: (Color online) Energy spectrum of the ground state 
(brown line) and first excited states forA'' = 50 molecules. The 
7i+ subspace is represented in violet and the Tii one in blue. 
All states in the same subspace are equally coupled to the 
field, and this coupling is different for each subspace, leading 
to energy levels crossings until the 7ii subspace becomes the 
ones with lower energy. In the inset, the crossing between the 
lowest energy state in subspace 7i+ (violet) and the lowest 
energy state in the Tii one (blue). 
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FIG. 2: (Color online) Logarithmic negativity Ld (blue, pink 
and brown dots) and L'^ (red and green dots) for different 
distances d between a pair of molecules and different positions 
p of a molecule in a chain, as a function of e^, for the lowest 
energy state with one rotational excitation. For blue dots, 
d= 1, for pink ones, d = 10 and for brown ones, d = 25. Red 
dots show Lp for p = 26 and green ones for p = 1. 



Cz for T varying from 0.2 to 1.2, which correspond to 
T varying from O.IK to 0.6K for KRb molecules. We 
observe that for all T, the entanglement decreases when 
ez increases. This happens because of the level crossings 
that exchange the energy ordering of the subspaces 7i+ 
and TLi. We recall that states in Tii are less entangled 
than the ones in the 7i_|_ subspace (see Figs. [T]and[2). 
We also notice that this entanglement field dependency 
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FIG. 3: (Color online) L26 a function of T and e^. 

is less pronounced with increasing T. This happens be- 
cause the population difference between each energy level 
becomes negligible, and we approach a situation where 
all entangled eigenstantes are equally populated and the 
system's state is separable. Finally, Fig. [3] shows that 
entanglement reaches a maximum as a function of T for 
all Cz- At this point, the excited levels are significantly 
populated when compared to the ground (non-entangled) 
state, and at the same time, there is a non-negligible vari- 
ation of population among them, so that entanglement 
properties are more pronounced. 

In the present paper we have focused on the first ro- 
tationally excited level, but our results can be gener- 
alised to highly excited subspaces. In such cases, sub- 
spaces with different entanglement properties and differ- 
ently coupled to the electrostatic filed can always be de- 
fined. Consequently, level crossings occur, modifying the 
entanglement and the physical properties of the system. 

In conclusion, we have made the first theoretical study 
of entanglement properties in a chain of polar molecules 
interacting by dipole force in the presence of an applied 
static field and with varying temperature. We have iden- 
tified two subspaces that can be labelled by their (AJ^)^ 
expectation values. Our results show a number of inter- 
esting and novel features, as a discontinuity on entangle- 
ment for the first excited levels when the electric field's 
amplitude is varied, irrespectively of the dimension of the 
chain. We have also investigated how pairwise entangle- 
ment depends on the distance between molecules. We 
showed that it decreases with distance, but persists for 
relatively far apart molecules, even if the interaction has 
been considered between nearest neighbours only. More- 
over, we have calculated entanglement of one molecule 
at different sites to the rest of the chain, showing that 



it also presents a discontinuity when Cz is varied. A 
striking feature of this system is that entanglement is 
absent at the ground state but increases with tempera- 
ture. We have thus studied how this entanglement tem- 
perature and field dependency takes place. The present 
study is relevant for fundamental reasons, since it con- 
sists of a many body system exhibiting peculiar entangle- 
ment properties, and also for experimental applications. 
The field of cold molecules trapped in optical lattices has 
made astonishing progress in the last years, and the re- 
sults presented here can be useful for the preparation and 
characterisation of entangled molecular states even at fi- 
nite temperature. The present work opens the perspec- 
tive of controlling entanglement properties of molecular 
chains by adjusting external parameters as the tempera- 
ture and external applied fields. 
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